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Abstract We report the results of our numerical simulation of classical- 
dissipative dynamics of a charged particle subjected to a non-markovian stochas- 
tic forcing. We find that the system develops a steady-state orbital magnetic 
moment in the presence of a static magnetic field. Very significantly, the sign 
^ ' of the orbital magnetic moment turns out to be paramagnetic for our choice 

of parameters, varied over a wide range. This is shown specifically for the case 



e 



of classical dynamics driven by a Kubo-Anderson type non-markovian noise. 
Natural spatial boundary condition was imposed through (1) a soft (harmonic) 
confining potential, and (2) a hard potential, approximating a reflecting wall. 
There was no noticeable qualitative difference. What appears to be crucial 
to the orbital magnetic effect noticed here is the non-markovian property of 
the driving noise chosen. Experimental realization of this effect on the lab- 
7 1 ■ oratory scale, and its possible implications are briefly discussed. We would 

like to emphasize that the above steady-state classical orbital paramagnetic 
moment complements, rather than contradicts the Bohr- van Leeuwen (BvL) 
theorem on the absence of classical orbital diamagnetism in thermodynamic 
equilibrium. 



C3 

on 



O 

o 



Keywords Fluctuation phenomena, random processes, noise, and Brownian 
motion — Diamagnetism, paramagnetism, and superparamagnetism 

O 

fT) . 1 Introduction 

The Bohr-van Leeuwen (BvL) theorem on the absence of classical orbital dia- 
£NJ ' magnetism in theromodynamic equilibrium has been a surprise of theoretical 

physics [1-3]. This BvL null result is statistical-mechanically exact. It is, how- 
ever, counter-intuitive inasmuch as a charged particle orbiting classically under 
the Lorentz force exerted by an externally applied magnetic field is equivalent 
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to an amperean current loop in the interior of the sample, and hence to a 
non-zero orbital magnetic moment. Moreover, the orbital moment is expected 
to be diamagnctic following the Lenz's law. A physically appealing resolution 
was advanced by Bohr [2] suggesting an exact cancellation, on the average, 
of the diamagnetic orbital magnetic moment in the interior of the sample by 
the paramagnetic moment subtended by the particle skipping the boundary 
in the opposite (paramagnetic) sense — the edge current [3]. More recently 
[4] , the role of the boundary was treated explicitly through a solution of the 
Fokker-Planck equation associated with the classical Langevin dynamics of 
the charged particle in a finite but unbounded space, namely the surface of 
a sphere (recall that, strictly speaking, a boundary has no boundary). Again, 
the orbital magnetic moment turned out to be zero. In our recent work [5,6] we 
have, therefore, re-examined the BvL null result. Our analyses strongly suggest 
that the vanishing of the classical orbital diamagnetism is a direct consequence 
of detailed-balance (the microscopic reversibility), namely that there are no 
cycles in thermodynamic equilibrium. Now, the detailed balance is, of course, 
conditioned mathematically by the second fluctuation-dissipation (II-FD) the- 
orem of Kubo [7] . In terms of the classical Langevin equation governing the 
stochastic dissipative dynamics that underlies equilibrium statistical mechan- 
ics, the F-D relation requires the driving noise to be markovian. Indeed, it 
has been shown [8] that a markovian non-equilibrium steady state can always 
be transformed to an equivalent thermodynamic state in equilibrium. It is 
our conjecture, therefore, that a stochastic dissipative system driven by a non- 
markovian noise may have in general a non-equilibrium steady state with hnitc 
orbital magnetic moment — without conflicting with the BvL theorem for ther- 
modynamic equilibrium. In what follows, we have addressed this issue through 
a detailed numerical simulation of the relevant stochastic dissipative dynam- 
ics involved. The results of our simulation support our conjecture, namely 
that the orbital magnetic moment is indeed finite for the non-markovian case. 
There is, however, a new surprise now, namely that the sign of the magnetic 
moment turns out to be paramagnetic! Morcvcr, for a classical gas of such 
charged particles with paramagnetic orbital moment, the inherently positive 
feedback may lead to an enhanced magnetic susceptibility — possibly even to 
a spontaneous ordering of the classical orbital magnetic moments. 



2 Stochastic dissipative dynamics in a magnetic field 

Consider the stochastic-dissipative classical dynamics of a charged particle 
(charge = — e, e > 0) and mass m in the xy-plane in the presence of a uni- 
form magnetic field of magnitude B directed along the positive z-axis. Let 
the particle be confined harmonically in the xy-plane. (In this simple model, 
the motion along the z-axis factors out). Harder confinement (reflecting wall) 
will be introduced later. The equation of motion can now be written down 
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F(t,At) 




Fig. 1 Shows schematically one possible realization of such a non-markovian noise compris- 
ing equi-spaced rectangular pulses of random gaussian height modelled on a Kubo-Anderson 
process. 



straightforwardly in the polar co-ordinate system as 

eB ■ 

mr = mrO 2 - kr - Tf - —rQ + F r (t,At) (1) 

c 

eB 

mr6 = -2mf6 - Tr6 + — r + F g (t,At), (2) 

c 

where, F r (t,At) and Fg(t,Ai) are the Kubo- Anderson-type noise terms [9]. 
The K-A noise, as depicted schematically in Fig.l, is best described as a 
sequence of random rectangular pulses of equal pulse- width At. The pulse 
heights are, however, identically, independently distributed gaussian random 
variates, of mean zero and finite variance. 

It is to be emphasized that At here is not necessarily a small quantity in 
any sense — it is to be regarded as a physical input parameter that makes the 
noise tunably non-markovian. The numerical simulation, however, must cover 
these not necessarily small time intervals At with much finer sub-divisions 
8t << At so as to ensure high numerical accuracy in treating the systematic 
parts of the driving forces present in Eqs. (1) and (2). More specifically, in 
our numerical simulation the number of finer time-step subdivisions has been 
taken to be > 500 (Figures 1-4). Moreover, the time of integration has been 
kept long (~5 x 10 8 sub-steps) in all these cases. 

It is convenient to define here the parameters J?o = confining harmonic 
frequency = yjk/m, fi c = the cyclotron frequency eB/mc, and 7 = r/m, 
the Stokes friction. Also, we write that F r (t,At) = mj 2 af r (t, At), and sim- 
ilary for F$(t, At). Note that a has the dimension of length, while f r ,fe are 
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dimensionless. Further, we introduce the dimensionless quantities ujq = ^o/j, 
uj c = Q c /l, T — it, and R = r/a the dimensionless radial coordinate. With 
these re-definitions, the equation (1) and (2) can now be re- written in the fully 
dimensionless form as 

R = R9 2 — LdQ 2 R — R — uj c R8 + f r (r, At) (3) 

R9 = -2R6 - R9 + u c R + fe(T,Ar), (4) 

where F r (t, At) = / r (r, At), and F$(t, At) = fg(T, At). Here the overhead dot 
denotes derivative with respect to the dimensionless time r. Also, / r (r, At) 
and Jo(t, At) are uncorrelated gaussian variates with mean zero and variance 
unity. 

The steady-state orbital magnetic moment M can now be written as 

M = {{-±(r x r)» = - (^pj ((R 2 9)) (5) 
giving the dimensionless orbital magnetic moment 




Here the double angle bracket denotes averaging over an ensemble of real- 
izations of the random noise as well as over time in the long-time limit. We 
identify fi with the steady-state value of the dimensionless orbital magnetic 
moment. 

In figure 2 we have plotted the dimensionless magnetic moment fi against 
the pulse-width At, for some chosen values of the magnetic field, measured in 
terms of the dimensionless cyclotron frequency co c . The parameters uiq and a 
have been set equal to unity. 

In figure 3, we have repeated, for the sake of comparison, the plot in figure 
2, but now for the case of a hard confinement. The latter is realized by modi- 
fying the soft, harmonic restoring force term — kr in the radial equation (1) to 
—kr(r/a) n , and choosing a large integer exponent n. The length parameter a 
here denotes the radial position of the effectively reflecting wall. With this, the 
term —u n 2 R in equation (3) becomes — uj a 2 R{R/A) n , with A = a/a. In figure 
3 we have chosen n=20, and A = 1. (Note that the limit n = corresponds to 
the soft (harmonic) confining potential). 

Finally, for completeness, in Fig. 4 we have repeated the plot (of against 
At) for the case of no confining potential, i.e., for oj — 0, with the other 
parameters remaining the same. 



3 DISCUSSION 

Through our numerical simulation of a classical model for the stochastic- 
dissipative dynamics of a charged particle, moving in a magnetic field and 
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Fig. 2 Plot of dimcnsionlcss steady-state orbital magnetic moment /i against the pulse- 
width At for some chosen values of the magnetic field (measured in terms of ui c ). This is 
for the case of harmonic (soft) confinement. For details, see the text. 
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Fig. 3 Plot of dimcnsionlcss steady-state orbital magnetic moment /i against the pulse- 
width At for some chosen values of the magnetic field (measured in terms of u} c = eB/mcy). 
This is for the case of hard confinement of nearly reflecting boundary as described in the 
text. 
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Fig. 4 Plot of dimcnsionlcss steady-state orbital magnetic moment /i against the pulse- 
width At for the same chosen values of parameters as in Figs. 2 and 3, except now for 
wq = (no confinement). 

driven by a non-markovian noise, we have demonstrated the appearance of an 
induced orbital magnetic moment in the steady state. Most significantly, the 
orbital moment turns out to be paramagnetid The appearance of a non-zero 
classical orbital magnetic moment in the non-equilibrium steady state implies 
deviation of the dynamical system from the second fluctuation-dissipation (II- 
FD) relation. The latter would have otherwise enforced a detailed-balance 
(meaning no cycles), and hence no induced orbital magnetic moment. The 
moment, of course, vanishes in the limit of delta-correlated white noise {i.e., 
At — > 0,(7 2 — > oo, with the product Ato 2 finite) that renders the system 
II-FD theorem compliant. The essential point to be emphasized here is the 
non-markovian nature of the driving noise used in the stochastic model. This 
point is consistent with the proven result of Prost et al. [9] that for a marko- 
vian dynamics, the non-equilibrium steady state can always be transformed 
into an effectively equilibrium state. It may be noted that our results hold for 
a soft (harmonic) potential confinement as well as for a hard (reflecting wall 
type) potential confinement. Also, the magnitude of the orbital paramagnetic 
moment is non-monotonic in the externally applied magnetic-field strength — 
initially increasing with increasing magnetic field strength and then decreasing 
in the high-field regime. Further, the magnetic moment scales as the variance 
of the random pulse-height for a given pulse width in the case of harmonic 
confinement (see Eq. 5). In the case of the hard confinement, however, there 
is no simple scaling, but the qualitative behaviour remains the same. 

For the sake of completeness, we also carried out a similar simulation for 
the case of no potential-confinement, i.e., with ojq = 0. Interestingly, but not 
surprisingly though, we obtained as in figure 4 a steady-state orbital magnetic 
moment which has the opposite sign — it is diamagnetic! This is, however, con- 
sistent with the physical picture of Bohr [1-3]: without the confinement, there 
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are naturally no orbits skipping the boundary, i.e., no edge current, which 
would have sub-tended a paramagnetic moment leading to the cancellation. 
We are thus left only with the Maxwell cycles well within the interior that 
constitute the amperean current loops giving the diamagnetic moment (the 
Lenz's Law). This avoided cancellation was seen also in the analytical solu- 
tion for the simple case of a markovian (delta-correlated) noise [10]. This now 
turns out to be true even for the non-markovian case as shown in our simu- 
lation. The overall picture is one in which the particle, initially at the origin, 
say, diffuses outwards on the average, but the orbital diamagnetic moment 
reaches its steady-state value on a relatively short time scale which is deter- 
mined by the parameters appearing in equations (1,2). The outward diffusion 
simply spreads out the total orbital magnetic moment (oc r x f ) over the ever 
increasing area covered, but without changing its time-averaged value — the 
steady-state diamagnetic moment. 

4 CONCLUDING REMARKS 

The paramagnetic sign of the induced orbital magnetic moment can have inter- 
esting consequences of considerable physical significance. After all, the para- 
magnetic sign of the moment inherently signifies a positive feed-back effect 
when we consider not just one but a system of many charged particles. Here, 
the mean self-field can in principle lead to a spontaneous macroscopic orbital 
magnetic moment. 

As for a possible experimental realization of such a confined system, we 
begin by noting that what is really essential for obtaining the classical orbital 
paramagnetism in a non-equilibrium steady state is the non-markovian nature 
of the stochastic forcing. Thus, a micron-sized sample of a semi-metal (such 
as bismuth) trapped in an optical tweezer, and irradiated with random laser 
pulses in the presence of an external magnetic field would constitute a possible 
candidate system. The laser impulses should impart high enough kinetic energy 
(high nominal temperature) so as to create a non-degenerate (classical) gas of 
charged particles (electrons and holes). Then, the high temperature washes out 
the quantum signature, i.e., the discrete quantum level-spacings in the micron- 
sized sample, leaving behind a classical charged particle system. (Recall that 
the orbital magnetic moment does not depend on the sign of the charge on 
the particle, e.g., be it an electron or a hole). It is then reasonable to expect 
the total induced orbital magnetic moment to scale up with the number of 
charged particles in the confined system. This mechanism may even generate 
spontaneously a macroscopic magnetic meanficld — such as the seed field of 
interest in the astrophysical context. 
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